function [s1,s2,s3,s4] = randdiv2(s,i,k,m)

    ss = s((i-1)*(3*k+m)+1:i*(3*k+m),1);
    s1 = ss(1:k,1);              % Size k for T at t=0  
    s2 = ss(k+1:k+m,1)';         % Size m for F1 at t=0
    s2 = s2(ones(1,k),:);        % Repeated k times for a kxm vector
    s3 = ss(k+m+1:2*k+m,1);      % Size k for F0 at t=0  
    s4 = ss(2*k+m+1:3*k+m,1);    % Size k for F1 at t=1

end
